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ABSTRACT 


The purpose of this thesis is to research the 
availability of targets moving on or near the earth’s 
surface when viewed by an orbiting satellite. A discussion 
of basic orbital mechanics is presented as well as a 
development of a suitable coordinate system. An analysis of 
non-linear observability is then provided. Lastiv an 


observer is designed and sucessfully simulated. 
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I. INTRODUCTION 


This thesis deals with satellite tracking of relatively 
low altitude targets. The long term goal of this project is 
to determine an effective general-coverage satellite orbital 
pattern. The short term goal of this thesis is to explore 
the possible use of various coordinate systems and to design 
a satellite observer. A target such as a slow moving 
aircraft is assumed for the basic development. 

Chapter 2 contains an overview of basic satellite 
mechanics and dynamics as well as a brief discussion of 
satellite detection and tracking urea Equations 
regarding orbital motion will also be introduced. Also 
chapter 2 provides an overview of the entire thesis. 

The theory and mathematics of satellite and target 
motion for a specific case are developed in Chapter 3. The 
geometry of motion on a spherical surface is detailed in 
that chapter. Specifics such as geometric swath width, 
line-of-sight and precession are furnished. This leads into 
an in-deapth examination of coordinate systems. Various 
attempts to find a suitable set of reference axes for 
analytic studies of the tracking problem are presented. Le 
is desired to achieve a target and observation model that 


has a convenient basis for analysis such as either linear or 


bilinear in state space. Several trial coordinate systems 
were developed mathematically but found unsuitable. One 
trial system is considered satisfactory under certain 
assumptions and is, therefore, completely developed. 

Chapter 4 presents a discourse on observability in 
general, and non-linear observability in particular. The 
concept of connectedness and univalence is discussed. The 
chosen system is then analyzed for observability. 

The topic of observers is presented in Chapter 5. A 
background on basic obsevers 1S provided. Finally an 
observer 1s designed and simulated by computer for the 
prefered system. 

Chapter 6 offers conclusions and recommendations for 


further studies. 


Pee tne oa le LL TE IN ORBIT 


Before an understanding of satellite tracking can be 
achieved it is essential to have an understanding of basic 
satellite dynamics and satellite sensors. The concepts and 
equations that govern the motion of a satellite in orbit 
around the earth are based on the physical laws put forth by 
kepler and Newton. These laws of motion apply to artificial 
satellites as well as planets and moons. The physical 
geometries and forces are the same. The following 
explanations and derivations form a basis for orbital motion 
of a satellite. 

A. THE SHAPE OF THE EARTH 

In reality the earth is neither exactly round nor 
exactly a rigid body. The constituents can be molded and 
deformed to a certain degree. The earth’s rotation causes 
it’s own materials to feel a force that pulls tangentially 
into space. Fortunately the earth’s materials also undergo 
a constant centripetal (center seeking) acceleration to keep 
them in place. This centripetal acceleration, which is 
directed towards the center of curvature, is provided by 
earth’s gravity. 

The gravitational force produced by the earth is 
directed towards the center of the earth. However, the 


momentum produced by the earth’s rotation pulls at the 


earth’s materials. Each particle of earth has both these 
forces acting on it. The earth’s surface at the equator 1s 
actually moving faster around the rotation axis than the 
surface material at or nearer the poles. Figure 2.1 
illustrates this. Therefore the materials on the surface at 
the equator feel this tangential pull the most. Analysis 
shows that the result is an earth shaped like an oblate 
ellipsoid. 

An oblate ellipsoid can be described as a sphere which 
has been compressed along the polar axis and therefore 
bulges at the equator. A cross sectional slice through the 
poles yields an ellipse. A slice along the equatorial plane 
yields a circle. 

This oblateness (measure of the earth’s flattening) of 
the earth is only slight. The equatorial diameter is 
calculated to be 12,757 km whereas the polar diameter is 
calculated to be 12,714 km. This is only a 43 km difference 
which is about 1 part in 297. 

B. GREAT AND SMALL CIRCLES 

The earth’s equatorial bulge can be temporarily 
disregarded in order to consider the earth as a sphere 
turning beneath an orbiting satellite. The intersection of 
a plane passed through the center of the sphere and the 
sphere is the largest circle that can be drawn on the 
sphere’s surface. This is known as a great circle. The 


shortest distance between any two points on the surface of 
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the earth is an are of a great circle. The earth’s equator 
is an example of a great circle. There are an infinite 
number of approximate great circles on the earth’s surface 
Since there are an infinite number of planes that can be 
passed through the earth’s center. 

There is a series of 24 special great circles that pass 
through the earth’s poles and are evenly spaced from each 
other. These are called meridians and they intersect with 
the equator at right angles. Meridians join at both poles. 
Forming other right angles to the Meridians are parallels. 
Parallels are small circles created by passing planes 
through the earth parallel to the equatorial great circle. 
See Figure 2.2. 

In order to determine precise locations on the earth’s 
surface, latitudes and longitudes are utilized. The 
longitude of a specific place refers to the are (measured in 
degrees) of a parallel between that place and the prime 
meridian (which passes through Greenwich, England for 
reasons of history). Longitudes run east and west. 
Latitudes, however, run north and south. The latitude of a 
apecific place may be defined aa the arc (in degrees) of a 
meridian between that place and the equator. Figure 2.3 
shows an example. When the earth is considered as having an 
ellipse as a cross section (which it has) instead of a 
circular one, the length of a degree of latitude is slightly 


greater at the poles than at the equator. 
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C. ORBITAL MOTION 

A presentation concerning orbital motion is basic to 
satellite tracking. Gravitational attraction and momentun 
combine to keep a satellite in perpetual orbit above the 
earth (Cin the absence of friction and external celestial 
attractions). 

1. Gravity 

The understanding of gravity and orbital motion can 
be traced back to Johannes Kepler (1571-1630) and Issac 
Newton (1642-1727). Kepler is known predominantly for his 
three laws of planetary motion. The laws are as follows; 
(1) Each planet moves about the sun in an orbit that is an 
ellipse, with the sun at one of the foci of the ellipse. 
(2) The straight line joining a planet and the sun Sweeps 
out equal areas in space in equal intervals of time. 

(3) The squares of the sidereal periods of the planets are 
in direct proportion to the cubes of the semimajor axes of 
their orbits. 

The above three laws apply equally well for 
artificial satellites as they do for planets. Newton 
restated and clarified Kepler’s laws. Newton was an 
advocate of rigorous proofs whereas Kepler preferred to 
State empirical laws based on observations. 


Kepler’s first law stated mathematically is 


f = mv7?/r C2 
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where: f=centripetal force needed for circular orbit, 


m=mass of planet (Cor satellite) in orbit, 


v = velocity of orbiting body, 
r = distance between sun and planet (or satellite and 
earth), 


Kepler’s third law restated mathematically is 

(m,+ mJp? = 4n?a°/G (2527 
where ml ander? refer to the masses of two bodies that 
revolve mutually about seach other. 
p = period of revolution, 
a = semimajyor axis of relative orbit, 


G = universal gravity constant (¢ 6.67 x 10 ka -5~ ) 


Newton also explained the reasons behind kepler’s 
original observations. That is, Newton supplied the laws of 
motion which are at the root of classical mechanics. 

These laws are as follows; 


(1) Every object remains at rest or in uniform motion 
unless an external force acts upon it. 


(2) The product of the mass of an object and its 
acceleration vary directly as the resultant force, and 
the change in motion takes place in the direction of 
that force. 


(3) For every action there is an equal (in magnitude) and 
opposite Cin direction) reaction. 


In the case of orbiting planets or satellites, 
gravity is the unseen force referred to in Newton’s first 
law. Newton also postulated that every particle of matter 


in the universe attracts every other particle. The force of 
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this attraction is proportional to the product of their 
masses and inversely proportional to the square of the 
distance which separates them. 
In equation form this is 

F = Gm, m, /d? G22 3) 
where d is the distance between the center of mass 1 and 
mass 2d. 

Orbital motion can be understood in terms of 
Newton’s laws previously stated. The critical factors in 
putting a satellite in orbit (by the traditional launch 
method) are speed and direction of movement at burnout. 
Burnout is when the rocket engine shuts off and the 
satellite behaves as an astronomical object. 

As a satellite follows it’s orbital path it 
continuously falls toward the earth due to the earth’s 
gravitational pull. However, the satellite’s momentum 
prevents it from really being pulled any nearer to the 
earth. The satellite’s orbit is the result of two main 
forces. The momentum of the satellite is a measure of its 
state of motion. The inertia of the satellite (recall 
Newtons first law) is that property that causes the 
satellite to resist acceleration and travel ina straight 
line. For the satellite to move in a circular path rather 
than in a straight line, it must continually suffer an 
acceleration toward the center of the circle. This 


acceleration is centripetal acceleration. The central force 
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that produces the centripetal acceleration is that 
gravitational attraction between the satellite and the 
earth. Figure 2.4 illustrates this in vector form. 

The period of an artificial satellite is directly 
related to the size of it’s orbit; 1.@, the closer the 
satellite is to the earth the faster it travels. 

2. The Center of Mass 

The center of mass can be defined as that point 
within a system that either remains fixed or moves as if the 
entire mass of the system were concentrated at that point. 
For a number of particles the center of mass is defined as 


Tiral 


: T ona C1/Mysrat 2.2m: F. (2.4) 


An example of this for a two body system is 


Kom = MX, + MX (27 
m+ ™m, 

Mass 11s a distance x(1) from an arbitrary origin, 
and mass 2 is a distance x(2) from the origin. XCcm) 
represents the location of the center of mass. A solid 
object can be thought of as a collection of a great many 
particles. In this case equation 2.4 applies. 

Very often it is useful in certain calculations to 
treat a particular body as if its entire mass is 
concentrated at one point. That one point is the center of 


Mass. 
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It can be shown (Ref. 1] that a large sphere 
attracts other bodies as though the sphere’s mass is all 
located at the center. This holds true as long as the 
sphere is of uniform density or is made up of concentric 
shells each being of uniform density. The earth is often 
approximated as being a sphere made up of uniformly dense 
shells (Ref. 2]. 


eye Equations of Motion 


A satellite in orbit around the earth will feel by 
far the strongest gravitational pull from the earth itselé., 
Other astronomical bodies, such as the moon or sun, do exert 
force disturbances but are too far away to exert a very 
strong pull on the satellite. Assuming the gravitational 
acceleration of the satellite is due only to the earth’s 
gravitational attraction, g can be defined as 

g = Gm/r? (2360 

The value of g does vary with respect to altitude 
and latitude (Since the earth is not really an exact 
sphere). Tables 1 and 2 depict the variations of g. The 
values of g are close enough to allow an approximation of 
the earth’s shape as spherical for most applications. 

An orbiting satellite is considered to be in the 
earth’s gravitational field. The vector symbol g refers to 


this field and is defined as 
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g = F/m C27 
When an object is a distance h above the earth’s 
surface it has potential energy defined as 
U = -C-mgd)h = mgh C225) 
Gravity is a conservative force pointing to the 
earth’s center and has the value (-mg). The gravitational 
force exerted on an object in this field can be derived from 


the potential energy equation as 


B=" au = 2 (-Sua] = -GMm (22) 
dr dr Te sag 


A more general form for potential energy invoking 
the universal constant G is 
U = -GMm/r a OP 
Potential energy can be converted into kenetic 
energy (CRef. 3]. As an object falls to earth it loses its 
potential energy as its height decreases but gains kenetic 
energy as its velocity increases. Kenetic energy is 


mathematically defined as 


K = GMm/2r = mvi?/2 C 2 isaie 
The total mechanical energy is 
E= K + J (252 
Orbital motion can be considered as motion in a 
plane. For one full revolution around the earth the 
satellite remains ina plane as it traces out an ellipse or 


circle. 


22 


The velocity, period and altitude of a satellite are 
all inter-related. The period of a satellite (Ref. 2] can 
be defined as 

C= 2w(R + ho/v , G2eaS) 
where v is the velocity of the satellite. 
The centripetal force, is often defined as 


F = mjv?/(R + h) (2.14) 


where m is the mass of the satellite. 

It is obvious to the casual observer that the 
centripetal force is also in part a function of the 
satellite’s altitude and velocity. Polar coordinates are 
used in figure 2.5 to illustrate the elliptical orbit 
variables. 

4. Uniform Circular Motion 

Sometimes a satellite’s elliptical orbit is that 
perfect ellipse the circle. If the satellite’s orbital 
velocity remains unchanged in a circular path then the 
satellite is moving with uniform circular motion. 

Since the orbit of a two body system Cearth and 
satellite) is often ina plane it is possible to use polar 
coordinates to develop basic equations. Figure 2.6 shows 
the relation between polar and rectangular coordinates. 
Velocity in polar coordinates for circular motion can be 


expressed as 


as 





Figure 2.5 ElL’ipeiecet Orbit 
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Figure 2.6 Polar and Rectangular Coordinates 
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Vv = ree (2.15) 
where 8 refers to motion in the theta (8) direction. 
Acceleration is found by taking the derivative of the 
velocity 

a = -r@?r + rae e216) 


Energy can be expressed as 


E=K + U = m(r*8*)/2 - GMm/r Geaa7> 
and angular momentum as 


L = mr’78 Ze 


Uniform circular motion is easier to deal with than 
non-uniform circular motion. There is only an angular 
component of velocity present in the equation. It is often 
useful to simplify the satellite dynamics in order to more 
easily manipulate equations. 

>- Equations of Motion in Polar Coordinates 

Consider the special case of a satellite ina 
circular orbit. This can be treated as motion in a plane. 
Polar coordinates are especially usefull for this case. The 
two variables that define motion and position are the radial 
coordinate, r, and the angular coordinate, 98, conversion 


between rectangular coordinates and polar coordinates is as 


follows. 
r= |X? + ¥? X = rcos(s) 
8 = arctan(cY/X) Y = rsin¢é8) 


Figure 2.7 illustrates this. 


2a 


Polar coordinates are often used any time 
curvilinear motion occurs. In Figure 2.8 basic curvilinear 
motion is shown. The differential length dr has components 
in both the r and 8 directions. The vector r can be 
expressed as rn where n iS a unit vector in the radial 
direction. 


Velocity can be defined as; 





v = d(r) = d(rn-) = dren, + rdn- (2a) 
- dt aie dt dt 





Rectangular coordinates are used once again to achieve an 


expression for the time derivative of the unit vector n-_ and 


Douce 

mn, = cos(s)i + sin(8)3 Ne = -sin(@)i + cos(8) 3 
n, = -Ssin(@)i8 + cos(8)38 = (-sin(9)i + cos(8)3)8 = Gn, 
No = ~cos(8@)i® - sint8)j38 = -(cos(8)i + sin(8)j3)8 = -8n . 


Velocity can now be expressed as; 
Vie rn. + rén , 
Following the same rules of differentiation acceleration can 


be derived. 


a=vV-=frn_+ r6n., + r8n, {SOT tO. nm 


—T -? 


Simplifying, 
a= (© - r@?7)n_ + (2r@ + r@)n, (2226p 
a tr 
For circular motion the radial component, r, remains 
constant. Only 86 varies. Therefore; 
v= r6n, (22a 


che ma 2 Os rén, (2.22) 
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The acceleration term directed in the negative radial 
direction is the center seeking acceleration. 
Figure 2.9 depicts the directions for the force vectors. 
The sum of forces is calculated as follows; 
SF = ma = m(r-r9?)n (2.239 
-_c ~—c —T 
SF. = ma_ = m(r@ + 2r8)n (2.24) 
— OQ ~eC 


Using differential equations and state form motion can be 


represented as in equations 2.27 and 2.28. 


dr/dt = V, (25255 

d@/dt = (2.269 

AVT7 at = Ve ee 7 ee or (25275) 

: aqea7dt = —Z2VeG/<5 + CUR Gy Za, (2.28) 
where U, = F./m and U,= Fo/m. Both variables represent 


thrusts. The symbols 6,- and 64 are used to represent 
disturbances which are suumed to be zero in this case. 
D. SATELEIGE Senses 

Various sensors on board the satellite are employed to 
keep the satellite itself on track, and to detect and track 
&@ possible target. The satellite guidance system uses 
sensors for measuring certain vehicle dynamic variables. 
Such variables include satellite acceleration, velocity, 
position and angular velocity. Guidance systems generally 
determine these variables, compare the received information 
with desired parameters, and generate correction commands. 
various types of sensors include inertial », optical and 


radio sensors. The satellite may make use of gyroscopes to 
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establish a set of reference axes on board. Active thrust 
control is used to align the body axes with the reference 
axes. Once the satellite has established the craft’s 
position and attitude active control components are 
activated. These components respond either to a telecommand 
from a ground station or to the satellite’s on-board 
computer in order to correct the satellites’s status. on 
some applications where perturbations are likely, feedback 
control is added to the active control of a satellite. In 
some cases optical sensors are used in conjunction with 
inertial Banh ite ree Optical sensors include such devices as 
sun finders, sun sensors, planet sensors, celestial 
trackers, and horizon scanners. The gyroscopes and 
accelerometers are examples of inertial sensors. In certain 
applications the inertial sensors provide short-term 
stability and optical sensors provide long-term stability. 
The gyroscope is an inertial guidance system that is 
often used for ships, aircraft and spacecraft. A gyroscope 
is a device that possesses a high rate of spin about an axis 
of symmetry that has freedom of angular rotation. Examples 
of gyroscopes are the rate gyro and the integrating gyro. 
Surveilance satellites normally have an entirely 
different set of sensors for target acquisition. Very often 
this type of satellite carries several types of information 
gathering devices. High resolution optical devices, radar 


scanners and infrared capacity can all be included in a 


Zo 


satellite payload. Passive systems are sensitive receivers 
which normally provide data on bearing and bearing rate. 
These systems have the advantage that they do not alert the 
target of their actions. Optical techniques normally 
involve lenses and cameras. This can be very precise, 
especially at low altitudes. 

Radar is an example of an active sensor. It bounces 
radio waves off a chosen target. This is excellent for 
determining both bearing and range. Both optical and radar 
methods can be degraded by poor atmospheric conditions. Fier 
the radar tracking technique noise is always the main 
limitation. 

Inverse scattering methods are being developed to enable 
polarimetric radar to obtain better target information. 

This procedure involves illuminating the target with 
polarized waves and observing the amplitudes and phases of a 
set scattered waves. 

Another type of sensor is the synthetic aperture radar 
system. This is based on holography methods. Holography 
has the ability to record three dimensional pictures and 
focus sharply on both the near field and the far field at 
the same time. This involves recording a wave interference 
pattern. A microwave generator is used to provide a 
constant frequency microwave signal as well as a reference 
wave. This is considered to be a highly precise and 


accurate sensor. 
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Passive sensors are normally infrared sensors. There 
are two major techniques used in enfrared detectors. The 
two main types are called photon detectors and thermal 
detectors. 

In photon detectors, the technique used involves a 
photon of infrared radiation being absorbed by a semi- 
conductor electron which raises its energy level into a 
conduction band. Photons with less energy than the band gap 
produce bo signal and are thereby effectively filtered out. 

Thermal detectors function be sensing the temperatures 
change resulting from absorption of infrared radiation by a 
suitable element. Usually, this absorbing element has some 
temperature sensitive electrical property such as 
resistivity so that the temperature change is sensed 
electrically. 

One device often used in passive detectors is the 
radiometer. This device is a broadband, dual frequency , 
low noise, solid state, remote controlled mechansim. 
Radiometers are used to obtain high resolution imagery from 
low altitude satellite’s. 

A promising type of detector is the Silicide Schottky 
diode based infrared camera. These are easy to manufacture 
and have excellent performance. The camera focal plane 
consists of a two dimensional array of metal electrodes 
fabricated on a silicon substrate. The focal plane is back 


illuminated. When an infrared picture is observed, hot 
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carriers are emitted from the focal plane electrodes. An 
electronix image of the scene is formed by accumulation of 
these carriers on a pixel by pixel basis. 

Another technique involves the Charge Injection Devices 
(CID). These devices are surface charge devices that 
collect photon generated charges and store them in MOS 
capacitors. 

Other types of infrared detectors include the Lead 
Sulfide (PhS) Detectors, Lead Selenide Detectors (Pb5e), 
Thermister Infrared Detectors and [Indium Antimode Charge 
Injection Devices. 

One tracking technique available to surveillance 
satellites involves the use of the doppler effect. The 
satellite emits a signal and then receives the signal’s 
return after it reflects off a given target. The change in 
frequency as a result of the relative motion between 
satellite and target leads to a calculation of the targets 
velocity. Ref 3 defines the return frequency, , detected 


by the satellite as 


v's Yl“ /7eo7 1 «uve, (Zacoe 
where 
Yos frequency detected if both were at rest, 
u. «= relative separation speeds, 


0 
iN 


speed of light, 
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Artificial satellites are often subject to a variety of 
disturbing forces. These include atmospheric drag, 
variation of atmospheric density, solar radiation pressure, 
surface charge drag, meteorite impacts, lunar or solar 
gravity caused perturbations and possible encounters with 
hostile killer satellites. These potential disturbances can 
act seriously to impede the position as well as the attitude 
stability of a satellite. Fortunately, there are passive 
and active stabilization systems to enable recovery from 


most perturbations. 
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III. SATELLITE COVERAGE AND COORDINATES SY > GEMo: 


This chapter includes satellite coverage, spherical 


geometry and coordinate systems. A coordinate system which 
enables analytical calculations of satellite observer and 
tracker is desired for at least a single case. Several 
systems are explored in an attempt to find a suitable set of 
reference axes to establish a base for more complicated 

A simple case is introduced and an 


analysis by computer. 


appropriate coordinate system is presented. 


A. GENERAL SATELLITE COVERAGE 
Basic information on satellite coverage is discussed 
here. An understanding of geometric swath width and 


precession is required for a detailed study of satellite 


tracking. 


(ee 





As a satellite travels around 
certain amount of surface area. This 
a spherical cap (figure 3.1). As the 
above the earth increases the area of 
increases. The area of the spherical 


calculated as 


A = 
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the earth it covers a 
area is in the form of 
satellite’s altitude 
the spherical cap 

can be 


cap (CRef.5] 


Cae 


where R is the radius of the earth and 98 is the angle shown 
a rigure 9S.2. 

It appears from this that a high altitude orbit is 
beneficial. However, it must be considered that as 
satellite altitude increases, satellite sensor accuracy 
decreases. 

ot fo target iS sat point “a” Cas in Figure 3.2) it 
can be viewed by the satellite as long as the satellite is 
between points 1 and 2 on its orbit. This portion of the 
orbit keeps the satellite above the horizon with respect to 
the target. The amount of time the satellite is above the 
Weoizon is. calculated (Ref. 5) to be 


t = 20/w C322) 


where w is the angular velocity of the satellite. 

As the satellite progresses in its orbit the spherical 
cap of coverage moves with it. As it moves. the cap traces 
out a ribbon around the earth. The width of the ribbon 
(Figure 3.3) is the geometric swath width (GSW). 

Figure 3.4 illustrates the GSW as the very dark arc on 
the earth’s surface. The dotted line from point 1 to 
point 2 is tangent to the earth at the surface location 
directly beneath the satellite. Ris the radius of the 
earth and h is the satellite altitude. The arc, a(l), in 
Figure 3.5 is one half the dark are in Figure 3.4. Simple 


geometry proves that arc atl) = RO. Therefore the 
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Figure 3.2 View of the Target 
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GSW = 29R. Figure 3.6 depicts 4 right triangle. rong. 
rules of trigonometry cosine(@) = R/CR+h). Therefore, 


8 = arccosine(R/(R+h)) 


It follows by substitution that 
GSW = 2RCarccosine(R/ (R+h) )) (3.28 
2. Ground Track and Coverage 
As a satellite travels in its orbit the distance it 
moves in one full trip around the earth is 
D = 2nr(R+h) (3.4) 
where R is the radius of the earth and h is the satellite’s 
height above the surface. The orbital speed is calculated 
(Ref.5S] to be ; 
So=4.2685(R/(R+h)) km/sec (Base 
The satellite also has a ground speed and a ground 
PEeack. The satellite’s position can be projected on to the 
earth’s surface. This point on the earth’s surface directly 
below the satellite is called the sub-satellite point (SSP). 
As the satellite moves its SSP forms a ground track which 
can appear as Simple as a circle about the earth or a mere 
point on the surface or they can appear very complicated. 
The ground speed due to the satellite’s motion alone is 
S = 7.9053¢(R/(R+h) ) km/sec (3.6) 
However, the earth does rotate at about 0.267472 
km/sec at the equator. The actual ground speed depends on 


the inclination of the orbit plane to the equatorial plane. 
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For all cases the ground speed is less than the actual 
orbital speed. As the ground track is formed so is the 
ribbon whose width is the GSW. Because of the earth 
rotation beneath the satellite the swaths may overlap at 
some places yet never cover other places. If the earth’s 
rotation is ignored and its shape is assumed spherical a 
satellite staying in one plane traces the same circular 
swath on every orbit. 

If the earth’s equatorial bulge is taken into 
account the plane of the orbit precesses about the earth’s 
axis of rotation. This precession is caused by the non 
central gravity force field in which the satellite travels. 
It ig measured by the rate of drift of longitude on the 
surface of the earth. 


For an elliptical orbit the rate of precession is 


Q = 9.95(R/a) « (cosine(id/(1-e?)? “GET 
where Q = orbital plane precession rate in degrees/day 
R = mean radius of the earth 
a = semimajor axis 
@e = eccentricity of ellipse 
i = angle of inclination of orbit 


For the special case of a circular orbit the 
precession rate is 
2 = 9.95(R/(R+h)) + cosine(i? (3.8) 
If the effect of this precession is ignored, the 


motion of a satellite lies ina plane passing through the 
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center of the earth. A satellite making one revolution per 
day at O° inclination to the equator is stationary over a 
point on the equator. In this case a stationary spherical 
cap on the earth’s surface is the only area covered. 
Commercial communications satellites are usually 
geostationary. This is a special case of the synchronous 
Orbit (24 hour orbit). If a synchronous orbit is inclined 
to the equator the ground track is a figure eight. 

3. The Spherical Triangle 

At this time all the simplifying assumptions are 

invoked to establish a base for further analysis. The 
satellite is traveling in a circular orbit above the 
equator. Obviously, transformations to other orbits are 
readily available. The SSP is always on the equator. There 
is a target traveling on a steady course at a constant 
velocity and therefore following a great circle heading. 
This is not an unreasonable assumption. It is normal for a 
ship or an aircraft to keep to a great circle heading at a 
most efficient cruising speed and altitude. The target is 
assumed to be on the earth’s surface or very near to the 
surface. Figure 3.7 illustrates the paths for both the 
surveilance satellite and the target. The target track is 
inclined i® to the satellite track. It is assumed that both 
can be considered to be moving on the surface of the same 


sphere with transformations readily available to the actual 
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satellite position. The earth is assumed to be a perfect 
sphere and is approximated by a point mass at the center. 
The two paths intersect in Figure 3.8 at point O. 

The letter 5S represents the present position of the 
satellite and the letter T represents the present position 
of the target. These three points (0,S and T) on the 
surface of the sphere define a spherical triangle. 
Spherical trigonometry is used to attempt to find a suitable 
coordinate system for the satellite tracking analysis. 
B. COORDINATE SYSTEMS 

Several coordinate systems are tested at this time. The 
set of reference axes is a relative one. It is assumed that 
the origin of the coordinate system 1S moving with the 
satellite relative to the earth. Appendix B explains. 

1. Use Variables $ and 6 

Supposedly any two independent variables from 

Figure 3.9 can form the basis for a state-space model for 
target tracking. The variables 2 and 6 are chosen for the 
first attempt to form a suitable coordinate system. The 
original four state variables x(1) through x(4) are defined 
below. Derivatives are found for the four state variables. 
In an attempt to establish a least-complicated structure, 
however, these lead to other state variables and other 
derivatives; 1.@.>, 


x€1)> = cosine(®S) 


x(2) = sine(Ss) 
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Figure 3.8 Intersecting Paths 
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x C3) 


cosine(d) 


x C4) s3ine¢s) 


Taking derivatives of two of the above variables yields 


x1) 


=-*% C2) 


x (3) -x (4) 


Using the trigonometric law of sines X(2) can be defined as 


x(2) =/sine(Msine(7)|= sine tt sine s) 
Sine (6) 1 Sine(a) 


Since in this case « and sine(a) are constants the second 
identity is easier to use. Sine ta) is replaced by the 
letter A and sine(&) is replaced by x(4). This yields 
; X(2) = sine(N)K(4)/A 
The same procedure is used for X(4). The law of sines 
provides an easy alternate form. X(4) can be defined as 
x(€4) = sine(T) A/sine( 9) 
Table 3 summarizes the mathematics of expanding the state 
space. It is soon obvious that this state space is not 
readily reducable to a single structure by increased 
dimenssion. It therefore is abandoned. 
2. Use of Variables 9 and $& 
Two new independent variables are now used in an 
attempt to generate a state-space model. As before the sine 
and cosine functions are employed. 


XC1>) = sinned) 
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TABLE 3. STATE SPACE FOR VARIABLES © and 6 


X€1)> = cosine?) X¥(€2) = sine) 
X¥(€3) = cosines) X¥¢4) = sinecd) 
XC1) = -X¢2) and xXC¢4) = -X¢3) 


X(2) = sin(f)sintS)/sinta) = sin()XC4)/A 
X(2) = cos(NNX(4)/A + sinM)X(4)/A 


Let X(5) = cos(TNX(4) and X(6) = sin(M)x<4) 


X(S) = -sin() (mp ?X(4) + cosMpnK(4) + cos(m MM) MX¢4) 


Let X(7) = sin(M)(M?7X(4) and = X(8) = cos(MpNKx (4) 


and X(9) = cos(N) (N) (NH) x¢4) 


X(7) = cos(N) ()”’ X(4) + sin) (2MNK(4) + sin(T CN) ?x¢4) 


Let X(10) = coa(n) (Ny X(4) and X(11) = sinh) (2MNXx¢4) 


and X(12) = sind) (MD 7KC4) 
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Ge > cosine($) 


X (3) sine(s8) 


X¢€4) cosine(8) 


The derivatives for two of the above are 


X C2) -x(1) 


K (4) 


-x (3) 

As before sine(a) is represented by the constant A. 
The law of sines is also used again. Table 4 summarizes the 
mathematics of the expanding state space. As before these 
variables as a basis for a coordinate system do not seem 
Suitable for tracking purposes. 

3. The Projected-Orbit Technique 

This technique involves projecting the circle which 
is the target track onto the equator which is the 
satellite’s ground track. Figure 3.10 illustrates the 
technique. The target position T is projected down onto the 
satellite’s orbit at position T’. 

Let + = ws) (the angular velocity of the satellite) 
and let 5 = w(t) (the angular velocity of the target). 

For convenience let Y = ST’ (see Figure 3.11) and X 
2 ST’ the arc length (see Figure 3.12). The arc length OT’ 


-N 


= OT and angle & = angle $’. It follows that 


6 = + + 8 
and 
8@=8- 1. 
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TABLE 4. STATE SPACE FOR VARIABLES & and @ 


XC1) = sine (s) X¢€2) = cosine(S) 
X(€3) = sine(8) X(€4) = cosine(8) 
X(2) = -X(C1) X(€4) = -XC3) 

KG ae sinony/-sinte) X(3) = sincar/ sinc@>) 

sin(a) Ssin(™| 


XC1l) = X(3)sin(D/A + cCos(TMTNXC3)/A 


7 Let X(¢5S) 


X(3)sin(M and X(6) = cos(mMnNx¢3) 


X(S) = X¢€3)sin@) + cos(Mnx¢3) 


Let X(7) = X(3)sin(N) and X(8) = cos(T)NX¢(3) 
X(7) = X(3)sin(Q) + XCB3rcostMN = X(Brcos(MN; (K(3) = 0) 
Let X(9) = X(3)cos(TN, 
X(9) = X(Z)sin(TY TM)? + X(3cosMN, 
Let X(10) = X(3)sin() (nN)? and X(11) = X(3)cos«APN 


° ae » 3% re 
XC LO) = ACSVeCos iy Gil Fees seem © 2 ry, 
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Figure S.11 The Line Segment ST’ 





pai 
Figure 3.12 The Arc ST’ 
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Taking time derivatives we find 
3 = z = x = W(t) - WCs) 

To obtain a state space interpretation Y and X and 
their derivatives are defined. Figure 3.13 shows an angle 
bisector which splits 98 and ST’. This forms two right 
triangles. Using geometry and trigonometry it follows that 

Sine(8/2) = T’A/R 
Al* =,S8772 
OT’ = 2Rsine(s/2) = Y 
Y = 2rcosine(8/2)8/2 = Rcosine(8/2)8 = 


Rcoosine (8-7) (wt) -wls) ) 


i X = ST’ = 8R ; 
X = OR + RO, 
but R = O; 


therefore, 
x = RE = RCwtt)-wls)) 

Recall R is a constant for this case and therefore R 
is zero. 

With Y and X defined above a logical state space is 
Originated. The measurement equation is also important. 

In order to define a measurement equation 
Figure 3.14 is constructed. Note that points O,S and T’ do 
not form a straight line. 

OT” 4 10S72 Si 


However, in terms of arc lengths the following applies; 


20 


—_—, 


OT’ = OS + ST’ 
The dotted line bisects angle «a and therefore 
BF = TT’/2 

Using geometry and trigonometry it can be shown that 

TT’ = 20OT’ sine(a/2) 
and 

BT/OT = sine(la/2) 

From Figure 3.10 it follows that | 

OT’ = 2Rsine( (g9+7)/2) 
Mena. triangle is defined in Figure 3.15. Using simple 
geometry it is obvious that 


OP ee AS ees 
and by simple trigonometry 


DT = TT’sineCf) 


and 


Dy TT’ cosine() 
tangent¢(€) = TD/SD 
therefore 
SD = TD/tangent¢€) 
The line segment ST’ is solved for 
ST’ = SD + DT’ = TD/tangent¢€) + TT’ cosine (fs) 
It is essential to express ST’ in terms of usable variables. 
ST’ = TT’ (sine(B)/tangent(€) + cosine) ) 
is achieved by substitution for segment length SD. 


TT’ = 20OT’sinela/s2) = 2¢2Rsinel(o+7) /2Z)sinela/s2) 


therefore 


a1 





Figure 3.14 The Line Segment ST 
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St aa eee geen Snes 2) sine (RB) + cosine(By 


aa _ tangent (€) 
Let 4Rsine((o+7T)/2)sinela/2) = Q@ for convenience. 


Then ST’ = Q@ ( sine (B)/tangent(é) = cosine(6)). Also, 
recall that 8 = WF/2 - a/2 


_—, 


ST’ = Q!sine(w/2-a/2) + cosine(nw/2-a/2) 
| tangent (é) es 


It is important to relate X and Y. 
X = ST’ = @R and Y = ST’ = 2Rsine(6/2) 
For convenience use go = X/R 
Y = 2Rsine(x/2R) and ‘en - X/o = Y/2Rsineto/2) 
therefore 
; X = oY/2sine(o/2) 
It is now necessary to solve for € which indicates bearing. 


Y = ST’ = Qcosinela/s2)cotangent¢€) + Asine(a/2) 


cotangent¢€) = iy - sinela/2) _ 


—— 


@cosine(a/2) cosine(«/2) 
Taking the arccotangent of each side, substituting back in 
for @ and then simplifying leads to 


€ = arccotangent. _sine(X2R) _ - tangent(a/2) 
 gsine(S/2)sinela) 


= 
This is the measurement equation. The final result is a two- 
variable state space and a measurement equation. From the 


measurement equation an observer can be modeled. Table 53 


Summarizes the results. 


a 


TABLE 3. @ STATE SPACE IN XAQAND S. 


State Space 
z = wit) (3223) 
X = Rewittrd-w d) C3. Tey 
Bearing information (measurement) 
€ = hieneaiial Sine(x/2R) = eg (Joli 
Sine(S/2)sine(a) 


Table & summarizes the results for an alternate 


coordinate system. 


TABLE 6. STATE SPACE IN Y AND &. 


State Space 
Ri- Y? Cw - w ) 


€@ = w 


Bearing information 
€ = arccet y - tan(ta/2) 
2Rsin(S/2) cosa) 





A previously researched thesis (Ref 6] provides a 
useful coordinate system based on latitudes and longitudes. 


Table 3.4 summarizes this coordinate system. 
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IV. OQOBSERVABILITY 


At this time it is important to determine observability 
of the target. Assuming the satellite’s detectors yield 
only the bearing information (Cas developed in Chapter 3), 
it is necessary to obtain the target’s position defined with 
X and 2. Since this cannot be measured directly it must be 
estimated in some way. Analyzing the state equations from 
chapter three it appears that if ~w+-and“,; are given the 
knowledge of X(O) and S(O) (the initialed states) are 
sufficient to determine X(t) and 3(t) all along the 
trajectory. 

X(t> = XCO) + Rf uyerr-uyrr3d 
° 
S¢Ct) = BCO) + ‘wr, (DAT 

2 
The possibility of estimating the state from the output 
measurement under the assumption that the system dynamics is 
completely known is defined as observability. Therefore 
the observability of the system must be determined. 

All changes in a given system must be reflected in the 
output for the system to be observable. Figure 4.1 shows a 
Simple example. The states X and x are estimated from 
knowledge of the output. If all states are observable for 
all times then the system is completely observable. This is 


defined in detail in (Ref. 6]. 


a2 


(X, ce) 
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Figure 4.1 A Simple System 


The test for observability on linear systems is fairly 
straight forward. In general the observability of linear 
systems is established by examining the linear 
independencies of the columns in the matrix function 
C(t St(t,t.), where C(t) is the measurement matrix and #(t,t) 
is the state transition matrix. For the linear time 
invariant system this leads to a simple rank test on the 
appropriate observability matrix. To illustrate this 
approach we may consider the following linear systen, 

K = AX 


Y 


CX 

where A 1S a square matrix, X 1S an n-dimensional state 
vector and Y is an m-dimensional measurement vector. 
Assuming that Yt) is differentiable up to the (n-1)th 


order, we have 


A matrix is set up as shown: 
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Y~ | = CA | x a yee 
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where Y is the extended measurement vector of dimension mn. 

The matrix @ is the observability matrix and has 
dimension mnxn. Clearly if rank @ = n, X can be uniquely 
solved in terms of the measurement Y. It should be noted 
that differentiability of Y(t) is not required for 
observability of linear systems and that the above result 
can be derived ina different way. The presented procedure, 
however offers the possibility of extending this result to 
the observability of non-linear systems; namely in the 


io 


equation, ¥ = Q@,X, any subset of n linearly independent 
equations can be considered Y = Q@X, @. can be regarded as the 
Jacobian J of Y with respect to X. Hence if rank J = n the 
linear time invarient system is observable. 
he TOR OERVABILITY OF NON-LINEAR SYSTEMS 

In the past it has been very difficult to determine 
observability of non-linear systems. Recently, however, a 
number of new methods for obsevability testing have been 
developed. 

A fairly simple test that determines observability in 
non-linear systems is developed in (Ref. 6]. This new 
method tests observability and identifies any unobservable 


states that may exist. If the system is unobservable than 


the measurement doesn’t provide enough information for state 


a7 


estimation. System observability and state estimation are 
directly related. 

A non-linear system can be represented as follows: 

X(t) = £(XC(t), ult),t) (4.1) 
where £() is an n-function, XER, UCR. 
The measurement equation is: 
Y(t) = NORCO), 2) (4.2) 
where h() is an m-function, YER. 

Y(t) must be differentiable up to the (n-1) order and 
UCct) must be differentiable up to the (n-2) order. A state 
Aa Cee is observable at t if knowledge of the input U(t) 
and the output Y(t) from time te to t , enables X (t.) to be 
determined. 

Equation 4.2 is differentiated to provide a necessary 


system of non-linear measurement equations. 





Y = h(x) 
Y=oh “On ox = ene CD 
st Fx Jt 
Y=2=jh+Jh, Jt + Shou = 
Oot ey ee Sula 


hie “oh ft Shee = eleGe up ) 
n-l 


Y = h {nz ) gb < hea) a < ey mee eee 


2 = t N-2 
cn ee 2) = h,., 6% ptt. Lec = -- Li ) 


+ 


Den-2) 


Define an mn measurement vector Y by 
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Y= = 
and an mn function H() to be 


H() = 
a 
The functional relation in vector form is 
Y = H(x,v?) (4.3) 
where v(t) is a function of yi PS meee er Cll 2) s 
With respect to equation 4.3 the question of observability 
can be developed as the existence of an inverse of this 
function. Clearly, this is related to the Jacobians of H. 
However, for the case of linear systems an inverse, if it 
exists, has to be unique. It will be shown late that here 
this need not be the case. 

Two conditions must be met for this non-linear system to 
be observable in the strict sense. These two conditions are 
connectedness and univalence. Connectivity is a necessary 
condition for observability in the strict sense anda 
necessary and sufficient condition for observabiiity in the 
wide sense. Connectedness is satisfied if every state is 
connected to the output in some way. Univalence is 
satisfied if every state is uniquely determined in terms of 


the measured output. A one-to-one mapping is considered 


a 


univalent. Both conditions together are necessary and 
sufficient to establish observability in the strict sense. 

To explain this consider Y(t) as expanded in a Taylor 
series. 

Y(t) = y(te) > y(t] (t-t > ery Ce «t-te 
* chy on-2) Poy Cte bet + r(t) 

Knowledge of the measurement trajectory Y(t) is 
equivalent to knowing the coefficients and remainder in this 
equation. This Taylor series expansion 1s considered unique 
and therefore the coefficients are also unique. The 
coefficients are the elements of the measurement vector Y. 
Therefore, any state X(t) is observable as long as it is 
connected in a one-to-one manner to an element of Y. if the 


connection is not one-to-one the system is observable only 


in the wide sense. An example of this follows. 


XS a 

In this example X can equal either the positive or 
negative square root of Y. This is a multiple valued 
function. Therefore the mapping is not one-to-one. Hence 
the system cannot be observable is the strict sense. 

As mentioned earlier the connectedness condition is 
satisfied if the existence of an inverse of the function 
can be established. According to (Ref. 6] the inverse 
function is considered established if the determinant of the 


nxn Jacobian of H does not equal zero for all x and v. If 
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J = 0 then one or more states are unobservable. If H has 
more than one right inverse then univalence is not satisfied 
and the system is only observable in the wide sense. 

However this will suffice for this case. 

B. JACOBIANS AND NON-LINEAR ANALYSIS 

Recall from chapter three that two possible bearing-only 
measurements have been found. These are illustrated in 
Tables 53 and 6. Both the measurement equations (bearing- 
only) are highly non-linear. The Jacobian matrices in 
Tables 7 and 8 are defined using the measurement equations. 

These Jacobians are necessary for observability 
analysis. It is essential to ascertain whether or not their 
determinants are equal to zero. The Jacobians in 
Tables 7 and 8 are based on the state space and measurement 
information available in Tables 5 and 6 respectively. 

The calculations for 7) » IDIS Ke NMG » SIS and J?)/\ 
concerning the variable X are in Tables 9 through 11. The 
result of the calculations for eee. S21 1g 8?) Vey y and 
INI S concerning the variable Y proved very complicated and 
are therefore not used. 

It is obvious that the determinants of these Jacobians 
are going to be difficult to analyze. Table 12 illustrates 
the determinant solution. To discover whether or not a zero 
value is achieved two separate methods are invoked. First, 
a Simplification is utilized. Recall from chapter one that 


a satellite is restricted regarding the amount of area it 
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can cover. Therefore the numerical value of the angles & 
and + are restricted. The variable X is composed of & and rT 


(X=R(S-7)). The development of the spherical triangle in 


TABLE 7 JACOBIAN UTILIZING X 





TABLE 8 JACOBIAN UTILIZING Y 
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TABEE 9 JACOBIAN - X - PART 1 
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TABLE 10 JACOBIAN - X- PART 2 
SH Lf x 7 
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TABLE’ Ti JACOBIAN - X - PART 3 
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TABLE La DETERMINANT OF JACOBIAN 
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chapter three assumes a limit to angles $&, +t and «. For the 
purposes of this work the range of possible values for these 


angles are confined as follows; 


$: 0 (plus epsilon) radians - w/4 radians 
tT; O (plus epsilon) radians - w/4 radians 
x: O (plus epsilon) radians - w/2 radians 


These limitations are reasonable for this particular 
case and they enable a useful Simplification. For small 
angles the sine of that angle can be approximated as the 
angle value itself Cin radians). Table 13 verifies the 
assumptions made here. 


TARDE St See SINE tr 1 e471 TONS 


a 


T/8 @/2 at $’s maximum 


r/8 -3926991 radians 


Sintnw/8) = .326834 





For this case the angle « is held constant. Using all 
the above assumptions and approximations the Jacobians and 
their determinants are greatly reduced in complexity. 

Tables 14 and 15 summarizes the development of the 
Simplified Jacobian and determinant for the measurement 
equation containing the variable X. This Jacobian is chosen 


because it is the least complicated of the two. 
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TABLE «14 


SIMPLIFIED VERSION 





| 
Measurements: | 


Te or S = sin? (a/2) 
Ra a = sin(a) 
7 


Jacobian: 


J =.¢ eee Ss Jeo 
‘ X? - 4RBSX + 4R757S Kix 
i =a, = ee ee 
: op X? - 4RSSX + 4R?57S~ J 


any XB - 8X 
(e - 4RBSX + 4R°B?S 
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| 

| 

| 

| 

| 
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— 
ul 


2X3 (X- 2ZRSS) + $(4R?57S- Kee 


\ U7 
Ah, * SRESK + GRE"? > aRSSK > 4R'E?S 7 





| 
| 


; fe 
J.2 = Te aR (X(4R?3?5S- XK?) +4RSXB(X- 2R)) 


Va DBR EU CD BGR cl (Wo =o acy linen = eee 
4RESX + 4R7S2S\ Xx? 4RESX + 4R?7S?75S / 
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TABLE 15 


LINEAR DEPENDENCE 





P= J J. = a?R? (XB - SX) 
X? - 4RE5X + 4R7S75 
since X? - 4RBSX + 4R*B*S = R*B%as 1 + on -~ bi? 
R¢a 
\ 
This holds for every X and & where a = sin(a) and 
Bem 25/4 = tanta/2>. 
\J{ = O if XB - SX = 0 (, - #£ 
“X £ 
tae Sees - tn © + Se nin) 2. = Aoe= SC 
‘Ss 
= @- = constant K =C,8 


‘ 


1 


| Therefore IS | 


O if X and & are linearly dependent. 
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—_ so elle Oe + a a 


o 


The second method involves the use of computer 
Simulation. The determinant (with the variable X as before) 
is calculated using no assumptions. The angles © and 7 are 
incremented in small steps from O radians to 1/4 radians. 
The angle aw (that is actually constant) has been given five 
different values for test purposes. The value of the 
determinant approaches zero but never actually achieves 
zero. 

Note from table 15 that, 


J | - Beta CAD = Sx) C4.4) 


et en ae = 


X?7 - 4RESX + 4R°S*s- 

If X@ = SX the determinant of the Jacobian is zero and 
the system is not observable. From Table 15 it is apparent 
that the system is not observable if X and $ are linearly 
dependent. 

The various calculations and computer simulations 
indicate that the system initially developed in chapter 
three (Table 35) is observable under general conditions. 
According to the results of the computer simulation the 
determinant never actually reached zero for the particular 
test values (50 test values) chosen. However there are 
certain conditions dependent upon the initial values & and 
7+, When the system ig not observable. Equation 4.1 defines 
these circumstances. Since the system is observable for 


most cases it is possible to design an adaptive observer. 


This 1s the goal of Chapter 5S. 
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V. OBSERVERS 


At this point in the satellite tracking project it is 
necessary to approach the state estimation problem. In 
general, a process can be characterized by the state and 
measurement equations; 

KC) = AX(t) Cont) 
Z(t) CxXc(b) Core 
where X(t) is an n-element column vector representing the 
states, d(t) is a q-element vector of measurements, A and C 
are nxn and qxn matrices respectively. The order of nis 
greater than or equal to the order of q. It is desired to 
obtain knowledge of the states of X(t). However, it i3 Z(t) 
that is measurable and therefore known. For this reason it 
is necessary to estimate X(t). The device for estimating 
the state X(t) is the observer. 
A. LINEAR LUENBERGER OBSERVERS 
A model of an open loop observer can be constructed as 


follows: 


al 


Y(t) = AYCt) G3:.3) 
This model operates as does the process or plant. The 
estimation error is 
Seen yrs X CLD (5.4) 
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This is differentiated with respect to time(t). 


VCE = ee 


e(t) = 
= AY(t) - AX(t) 
= AeCct) (3.20 
Therefore, 
a(t) =e e(O) = 8(te(O) (5.6) 


If the eigenvalues of A all have negative real parts 
then e¢t) approaches zero as t approaches infinity. The 
error signal ett) decays at a rate determined by the 
location of the eigenvalues of A. The rate of response of 
the process is also determined by the eigenvalues of A. 

To make the decay rate independent of the dynamic 
process, often a state model is ’driven’ by an error signal. 
In this case the observer is characterized by; 

Y(t) = AY(t) + GCZCT) = Gre) (S272 
and the error equation is 
eve). = YC) =Facee 

The matrix G is an nxq gain matrix which can be selected 
to determine the rate at which yt) approaches X(t). The 
observer equation can be written as; 

Y = AYCt) + GCCXCED=CY(t)) = AYCG) + GCeme (ScD 

The driving term is GCe(t) and its purpose is to drive 
the estimate Y towards X. The error equation is 
differentiated with respect to time. 

aCe Ce mee 


ect) 


AYCt) +*GCGKCt? -— €CxXCE)) — Axe” 


2 


eCEe- CA-+GORCY (tO CE) > CS o 


e(t) = Fec(t) (5.10) 
where F = A-GC 
e(t) = e - e(0) (S711) 


As expected, the results depend on the initial 
conditions. The gain matrix G must be chosen to place the 
eigenvalues of F at suitable locations. In designing an 
observer it is usual to place the eigenvalues of F where 
desired and then determine G. 

If the process has an input signal, u, then the 
appropriate equations are; 

; X(t) = AX(t) + Bue (S02) 

Z@e) ~= CXCED 

The observer for this system is characterized by the 
following equation. 

Y(t) = AY (t) +BU(t) +G( Z(t) -CY(t)) CaS) 

The error response is the same as in the above and is 
aiven by equation 3.11. A block diagram of the observer 1s 
shown in Figure 5.1. 

When designing a linear or non-linear observer it is 
desirable that the observer error become small rapidly and 
that the observer not be very responSive to noise. It is 
very difficult to meet both these goals since they seem to 


eont Li ct. The observer must also be supplied with a set of 


initial conditions. 
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If some of the states of a given system can be measured 
there is no reason to estimate them. An observer of reduced 
dimension can be designed to estimate only those states that 
cannot be directly measured. For example, if there are n 
states ina system and q of them can be meaSured, the 
observer only requires (n-q) states. Figure 35.2 illustrates 
the reduced order observer. 

B. NON-LINEAR OBSERVER 

For the design of the non-linear observer the techniques 
developed in (Ref. 7] can be applied. 

Declare two given n-dimensional vector valued functions 
of time to be Y(t) and 86(t). Let these two functions pe 
related by H(8) as Y(t) = H(8). The vector H(S8) is 
considered known. 

Assume that each ith row H(8) is differentiable at least 
once on all parameters 8. Therefore an nxn Jacobian matrix 


is defined as 








3¢o) =J HCO) = | 3 c2)| (Seep 
a J 
jo 8 
J;- (8) =JSHi (9) (5.155 
Also; 
Y(t) = Hea) + JHce) Je = Jca)act) «5.76 
5S 8 SH 
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Figure 5.1 Basic Observer 





Figure 5.2 Reduced Order Observer 
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Figure 3.3 is the general non-linear observer referred 
to in (Ref. 7]. The transpose of the Jacobian is represented 
as Toon A string of parallel integrators are denoted as( 
and 6,(t) represents a solution of Y(t) = H¢8) for any given 
trajectory ¥, Ct). In regards to this figure there exists a 
positive scalar s and a time T>O such that for time t>oT 

(l] 8s Ct) = 8, ¢t)>|| < s) if the following three conditions are 
met. 

1. K is positive definite 

2. The magnitude of the determinant of the Jacobian is 
bounded both from above and away from zero from 
below for the particular trajectory Y, ¢t). 

3. The magnitude of Y, ¢t) is bounded from above. 

In addition, by increasing the minimum eigenvalue of K 
the scalar s is made arbitrarily small. 

An estimation error signal is defined as; 


e(t) = 6;¢t) - 8,(t) (Seize 


Rearranging provides; 


ul 


Q@;¢t) St) 9+ TEE» 
Referring to Figure 3.3 the following relation is 
stated: 
@, = -KJ'(@,) (H(@,) - H(@,)) (5.18) 
For notational convenience the following is defined; 
=(e,8,) = Hte+8,) - H(@,) = HC8,) - HCA.) Cone 


In this case sigma (2) does not refer to a summation. 


It is simply a variable. 
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Figure 5.3 Non-Linear Observer 
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It follows from the above calculations that: 
e= 6, - @ = -KJ'(@,) =e,8,) - 4, (5.20) 
The time-varying Lyapunov functions of the error signal 
are defined as: 
Vie,t) = %F (e,8,)£(e,8,) (5.21) 
Since the determinant of the Jacobian is bounded 
2€e,8,>) = O if e = 0 them 12 follows (fnatle, (ee, =) om 
e = 0. Further calculations are summarized in Table 16. 
The minimum eigenvalue of i orar cer cno is denoted by 
Ne which is always greater than zero. The minimum 


eigenvalue of K is denoted by ae which is also greater than 


zero so both eigenvalues are positive. It follows that; 


Sealine 
v< - Ne \ECe, 8, 2-4 Ta, ECE, 84? (5.24) 
Y, = J3(@)) Oy (5.25) 
. =i ‘ 
@, = J(8,> Y (5.26) 


Since Y and J(8 ) are both bounded scdy and a, are 
also bounded for all t >» O. Therefore Ou < c, and 
[FCB] < So. 
Vi <-ApAcHECe,O 9? + CC. [!ECe, a, >} 
Vi< -CiNpde [ECe,8y 1-C.C27y 2A] 2740, C2/2A Ac (8.27) 
C. DEVELOPMENT OF SATELLITE OBSERVER 
For the design of the satellite observer (the system 
initially presented in Table 5S) the state 8 as described 
above now represents both X and $ as equation 3.28 shows: 
6 = {X,| = | X (3. 25> 


x 8 
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TABLE 16 LYAPUNOV FUNCTION OF THE ERROR SIGNAL 


Viit)= JV (et) Jom SV (et) 
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oe 
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The results developed earlier in this chapter 
(part B> are now applied to the particular problem of the 
satellite tracker. The Jacobian of Table 7 is redefined in 
equation Ss. 24.1. 

It 1S intuitively obvious that; 

Y = HCX,8,U) = HC8) 

The further development of both H(8) and 8 is presented 
in table 17. The information in Table 17 leads to the 
design of an observer for the satellite. This design is 
depicted in Figure 5.3. Figure 5.4 shows the model used to 
check observer results. Both the observer and model are 
Simulated using the FORTRAN programming language. The 


results are discussed in the next chapter. 
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Aber, 37 CALCULATIONS FOR OBSERVER “DESIGN 
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Figure 35.4 Satellite Observer 
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VI. SIMULATION RESULTS 


The basic observer designed in Chapter 5 is simulated to 
prove authenticity. Recall that the original state and 
measurement equations were developed in Chapter 3. Due to 
the satellite’a limited GSW, several parametera (including 
anglea 2 and 17) are confined in their maximum attainable 
values. Therefore, these parameters were kept within 
specific ranges. The angle alpha (a) is taken to be a 
constant 30°. This value is selected for convenience. 

The angular velocity of the target (*.) is chosen to be 
0.0005 rad/sec. This sets the target speed at about 
2mi/sec. The angular velocity (uU;) of the satellite is 
chosen to be 0.00045 rad.sec which yields a speed of about 
1.8 mi/sec. These are arbitrary but not unrealistic values. 
A. THE BASIC CONTINUOUS OBSERVER 

Figure 6.1 depicts the basic simulation block diagram. 
The initial conditions & and X are 0.2 radians and 1000km 
(621.5 mi> respectively. These values are arbitrary but 
within the acceptable range. Recall that X = R(&-7T). 

Figure 6.2 is a block diagram depicting the basic 
observer design that is simulated. Note that Figure 6.2 is 


virtually identical to Figure 5.5. 
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Figure 6.1 Rasic Simulation Bleck Diagram 





Figure 6.2 Observer Bleck Diagram 
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The equations representing the various elements in 
Figure 6.2 were initially developed in Chapters 3,4 and 5. 
Table 18 summarizes these equations. 

The gain matrix (as developed in Table 17) is chosen to 
be a diagonal matrix of the form: 

K = - 0 

Olas 

The elements k and k both must be greater than zero. 
A diagonal matrix was chosen because it is fairly easy to 
work with and it allows the error changes to be monitored. 

There are two ways to describe the observer to be 
Simulated. It can be described in the continuous domain or 
in the discrete domain. 

The observer was initially represented as a continuous 
model. A forth-order Runge-Kutta method is used for 
integration. The continuous model was used in order to 
eliminate any possible errors that could be introduced as a 
result of discretization. However, simulations for both 
versions produce no real differences. 

Be THE DISGREPE OBSERVER 

With the advent of micro-computers a discrete observer 
acting as a one-step predictor is a more realistic choice. 
In this model the foreward difference approximation is 
employed. 

The basic definition is; 

X(t) = X(t) 4 at KCK) 


rst] 
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VAG L Ewer OBSERVER EQUATIONS 
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1. Full Non-Linear Observer 

First the original, highly non-linear observer 
(figure 5.5) was simulated using a FORTRAN program on the 
IBM 370 system. All attempts to find values of gains and 
integration step for which observed states converge were 
unsuccessful. Typically the observer enters steady state 
with very small errors e, and e.3 however, the estimated 
states were far from actual values. It is possible that 
this behavior is the consequence of the sinusoidal periodic 
functions (sin(x/2R) and sin(t$/2)) in 7 and 7). Recall from 
Chapter 4 that this observer is not shown to be observable 
in the strict sense. Univalence may not exist due to the 
reiterative nature of the sine function. In order to avoid 
periodicity an approximation of a non-linear observer is 
considered. 

The simplified observer has the same form as the 
full observer (see Table 18), but the measurement equations 
and the corresponding Jacobian are approximated as developed 
in Chapter 4 (and repeated in Table 19). 

The simplified observer simulated was being very 
insensitive to changes in X and &. The observer stabilized 
with a small error (Ce, and es) while x and @ remained 
Significantly different from X and & respectively. The 
exact reason for this behavior is unknown. However it is 


theorized that the relative smoothness of the arccotangent 


curve Carccot 77) is responsible for this behavior. 
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TABLE (2 SIMPLIF LEDMOBSERVER 


Measurement equations: 


set s = sin? (a/2) and a= sin(a) 


h. Gigee> arccot [x ~ = 
RZa 


h,(X,B) 


wah (- XB = BX 


Jacobian: 


Jv = VOLIX = -aR 
X? - 4RBSX + 4R?575S 
Ji =) 7/SB = aR 
X? - 4RBSX + 4R?7575 
Tis, = SY/SX = aR 
X? - 4RBSX + 4R7575 
2XB(X-2RES) + B(4R?S?S-X?) 
X? - 4RSSX + 4R?737S 
Jy ee aR 


X(4R?B2S-K?) + 4ESXB(X-2RB) 
“Ya a aR Sie Smear care 
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\ 
X* - 4R@5X + 4R*875 / 


Soc 29 eee 
X? = 4RESX + 4R72¢5 


*® 


® 


4 
i 


| 


In order to circumvent this problem the measurement 
Signal is redefined. 


2. Redefined Non-Linear Model 





A new measurement variable is defined as: 

AY Ct) = x(t)/BCt) = Rsin(a) (cot ¢(7(t)? + tan(a/2)) 

If the value of «a is known and 7/(t) is measured; 
then uCt) can be calculated. Taking into account the 
assumed range for S(t) and X(t) and the values for 
parameters R and «a, the expression X/(RSsin(a))-tan(a/2) 
never achieves a very large value. Therefore “/ (t) is 
bounded away from zero. This implies that cotan(/(t)) has a 
bounded value. 

The derivative of p(t) is: 

L(t) = Rsinta) (= to/sin* (Z(t)? 
where 7(t) and 7<t) are measured values and £.<t) is 
evaluated. 

The redefined non-linear model and observer are 
presented in Figure 6.3. It is noted that this observer has 
the same form as the observer given on the block diagram in 
Figure 6.2. Observer equations are shown in Table 20. 

All attempts to find the power K,, K, and At (‘the 
time increment) that would force the states to converge 
failed. A satisfactory theoretical explanation for this 
behavior was not established. However, one possible 


explanation is that e ande are extremely sensitive to 
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TABLE 20 


MEASUREMENT EQUATIONS 


REDEFINED NON-LINEAR MODEL 


Measurement equations: 














h, (XS) = LOX pe) =. 
S 
h2(X,8) = -X,8) = XB - SX = 4 (X - SX 
52 S S 
Jacobian: 
Jy, = he = 1 re Ae | 
oS X 5 SB % 2 BS | 
; | 
Jz, = dhe = -10s) 
J X B? 
Jip = Dh, = XCB* - 2BCKB-BX)” = -1 (xB + 23x | | 
+ g 4 $3 : 
eS oy 
a(S 
Observability; 
[J = J, Jo ro Ji. Jz | <7 1 / 28x en i mi Bx = 4, /8x = *) 
S| & $? §& $e‘ 
J i> 0 1 f tye 2 x, 


x Ms 
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change in states which results in an extremely slow 
convergence with a highly oscillatory transient response. 


Ss Redefinition of State Variables 





Analysis of the model with redefined measurements 
reveals that the form of the equations allow redefinition of 
state variables. This may reduce the sensitivity of the 
model. 

A new state vector [U,V] is defined as: 

U = X/S and V = 1/8 
Table 21 illustrates further development. 

The block diagram in figure 6.4 shows the discrete 
observer. A simulation of the discrete observer was run 
with “k= 0.0005 and “$=0.00045. Several different gains 
were used to test convergence. Convergence is obtained for 
the following ranges: 

Oe< K,4t < 2 and O< K2At/|J22 | <2 


A ”~ 
Arbitrary initial values for U and V were employed. 


It seemed natural to assume Ue = O, X = 90 and Vo 1 (2 =1 
ia within the allowable value for ®). 

To examine the observer’s ability to track varying 
angular speeds it is assumed that Ujand “: are sine waves. 


rn 
CIT A sin(Dt) 


a 
U/s = B sin (Dt) 
Different amplitudes and periods are tested and 


successful tracking is achieved. Appendix C contains the 


resulting computer graphs. 
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Figure 6.4 Discrete Observer 
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TABLE 21 


OBSERVABILITY CONDITIONS AND DISCRETIZED OBSERVER 





Measurement equations: 
AMKU,V) = hy CU,V) = U 
ALCU,V) = h (U,V) = VCX - BUD 


Jacobian: 


JU OV 
Teen. GeCSVemee. 25h: = X - SU 


- 


Observability conditions are the same: 


ee - - | 
i = eee = 0 X- 8X #0 XA Bs | 
he | 


Continuous observer: 


“~~ 


U = K,J, @, + K,J., ©, 


* 


V 


K 2J,, [a K2J.2 @2 


SE TE el 


The discretized observer becomes; 
“ ine an 
— ba CUrc ec? — h, CUK eV.) Ho Ue 
: “ oN ‘ et ’ tae 
22,2 CU. ,V, ?D — fecu, » Vee Saye V_ CX - @ Uy) 


a ow “Aw 
U_,, = Un + K,atew + K,at s J2,Un,V Qe, , 


an 
Kel U 
7 “~ rx “Ww 
Vic ry = Vax 7 K,atJ22 (Un VOder2¢ 


——— es 


Js 


As expected the amplitude of the sine wave 
representing state (A/D;B/D) which can be tracked is 
inversely related to the period of oscillation. 

4. Model Decoupling 
Analysis of the error equations provide; 
e, = /ACU,V) - AKU,0) = U-U 
e. = LUU,V) - hiU,V) = V(X-8U)-VCX-8UD 
= X(V-V) -BCUV-UV> 

The block diagram in Figure 6.4 indicates the 
possibility of decoupling the observer into two first order 
models by letting J:, equal zero (Jiz2 is already zero). 

i Figure 6.5 shows a block diagram of the discrete 
subsystem. The discrete system of Figure 6.6 is certainly 
stable if it’s eigenvalue lies within the unit circle. 

l1-a} € iy 0<a< 23987 < K At < 2 

Under this condition the system has a steady-state 
error of zero for a constant input u. 

Hence the system acts as positioning servo-mechanism 
and is capable of tracking the relatively slow varying 
inputs with very small errors. 

3S. Subsystem Two 
Assuming that subsystem 1 is tracking u ¢ uct) = 
uC€t)) the error @,can be approximated as: 
; a iN 
e2.* (X - BSUICV - V) 
The observer is reduced to the gradient algorithm 


where the gradient is: 
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Sugcre ts. & Discretized Bleck Diagram 
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Figure 6.7 illustrates the observer. This result is 
also in accordance with the method proposed in (Ref. 6]. 

Defining V= (X-BIDV as a new measurement the 
corresponding function h is h=(X-BU)V. The Jacobian is 
ys SiufpSve x-8u. 

In order to simplify the stability analysis the 
Jacobian, J, is replaced with the signum function, 
Signum(J). This is sufficient to guarantee the negative 
feedback of subsystem 2. Subsystem 2 is analogous to 
subsystem 1 and will be stable if O< K, dt < 2. 
Consequently this system also performs as positioning servo- 
mechanism and will track relatively slow varying inputs V(t) 
with very small error. 

The presented analysis leads to a definition of the. 
decoupled observer. This observer is depicted in Table 22. 
Appendix D contains the computer simulation results. 

Figure 6.8 is essentially the same as Figure 6.4 
with Jzi = O and Jrzz = sSsignum(J22). Therefore it is 
possible to consider this algorithm a simplification of the 
algorithm developed in (Ref. 7]. However, it should be 
pointed out that the decoupled observer can be derived 
directly through analysis of the two linear subsystems 


presented above. 
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YABEE 22 


DISCRETIZATION 
ss wwe 
Seay U = - h,(U,v) 
‘ a 
= = _ = eee 
e. VCX SU) i h,(uU,V) 
U = Ke, 
4 
V = K,Sygn(J220¢e, 


~ 


After discretization this becomes 


f> “~ 
So ot ALS 1 
Vig, = Ve + K, at sygn(J2.20¢2. 
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In the decoupled model it is easier to control gains 
Since the J,., element does not affect stability. It is 
noted that slightly more oscillatory behavior is present 
at the beginning of the simulation. This is not unusual for 
gradient type algorithms. 

Satisfactory results for the product Kat = 0.1 are 
expected if the observer is analyzed from the sampling 
theorem point of view. This presupposes that both 
continuous subsystems posses only one time constant “= 1/K. 
The sampling theorem requires ane t < (/2 and atK < 0.5. 

Experience has shown that the best results for the 
one-step integration are obtained if the integration step 
(sampling time) is ten times less that the smaller time 
constant in the system. This means 4t * 0.1¢1/K) and 
itK = O.1. 

The sampling theorem is violated if 4tK=1. Asa 
consequence the observer requires much more time to extract 
sufficient information from the measurement and to start 
state tracking. The simulation results demonstrate that 
when the product AtK equals unity convergence starts after 
fifty seconds. This iS approximately ten times slower trom 
previous experiments where the sampling theorem is 
satisfied. 

In the case of 4tK = 0.01 the system pole is very 
close to the unit circle. Therefore, convergence is 


relatively slow. 
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As demonstrated by the various computer generated 
graphs, the simulation results for this system are 
satisfactory. Therefore the selected coordinate system is 


proven acceptable. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


The objective of this thesis is twofold. The first part 
of the thesis involves finding a suitable coordinate system 
in which to establish a satellite tracking model. The 
second part of the thesis uses this coordinate system to 
obtain an observability analysis and design an observer. 
This latter portion is the main thrust of the thesis. 

The satellite model is based on the coordinate system 
developed in chapter three. The model is simplified yet not 
unredlistic. The ma@in limiting aasumption confines poth the 
Satellite and target to a constant speed and heading. An 
unchanging aie oe ee is important due to the desire to 
ma&intain a constant angle a. These restrictions impased 
on the satellite are explained in detail in chapter three. 
Satellites and targets (aircraft) would normaliy travel a 
great circle path which requires a constant bearing. 
Therefore, the assumption 183 valid. This basic premise 
resulted in the epherical triangle that forma the bae2#1ag of 
the system presented in Table 3. Once the basic model 
became established an observability analysia was 
accomplished in chapter four. The basic observability 
theories developed in (Ref. 6] ‘(for mechanical springs) has 
been successfully applied to satellite observability. The 


system proved to be observable in the wide sense. 
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Further studies in the area of satellite tracking and 
observability are highly recommended. A more realistic 
satellite model could be devised in these studies. To 
improve the accuracy of the satellite motion equations the 
earth must no longer be assumed spherical and the 
Satellite’sa orbit must no longer be confined to simple 
Circular orSare-. This leads to a different coordinate system 
by necessity. Relative coordinate systems should be 
explored in more detail. Avpendix B discusses this Driefiyv. 
A possible approach involves using the Euier equations of 
motion. A target position with respect to a sateliite prime 
axils could be transtormed via pitch, yaw and roli to a 


es 


b 


satellite normal axis. Enhanced satellite tracking stud 
should ee a target that is not confined to a constant 
heading. The target should be allowed to alter course to 
further complicate the coordinate system and observer model. 
Both the satellite and target should be subject to random 
disturbances. This leads to the development of stochastic 
models. References 6 through 8 are highly recommended to 
anyone pursuing this course of researcn. 

This thesis has successfully developed a model for 
satellite tracking. The system has been determined to be 
observable and an observer has been designed. The observer 


model has been successfully simulated. 
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APPENDIX A 


EPENENTLS OF ORBIT 


Satellites travel in elliptical orbits. Figure A.1l 
depicts a basic ellipse. The distances a,b andc are 


related by equation A.l1. 


Cow 422 =>. CA) 
In this equation: @ = semi-major axis 
b = semi-minor axis 
c = distance between foci 


Eccentricity is also a parameter often used in 
connection with elliptical orbits. PCECHhEELICLUy, €, 15-4 


meassure of the elongation of the ellipse. Equations A.2 


and A.3 relate eccentricity to previously defined 


parameters. 
e? = 1 - (¢b/a) GA 2) 
c = ae (Agee? 
Other elements include apogee and perogee. Apogee 15S 


that point on the orbital ellipse where the satellite is 
furthest from the earth (focal point). The perigee is that 
point where the satellite passes closest to the earth. 

These new parameters are related to a and e by equations A.4 


and A.5. 


NN 
Oo 
A 

il 


apogee distance a(l+e) CA.4) 


perogee distance aci-e) C(A.5) 


I) 

QO. 
~y 
J 
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Any two of the above 31x parameters can define an 
ellipse. For the special case of a circle; e = 0, d =d 
c = O and a= b. 

To describe satellite motion more information is needed. 
The inclination angle, 1°, is the angle of intersection 
between the orbital plane and the earth’s equatorial plane. 
Figure A.2 illustrates. Intormation on the period of 
revolution of the satellite is needed as well as knowledge 
of the precise time that the satellite passes it’s apogee or 
perogee. 

Some references use information on the ascending node to 
determine satellite paths and positions. Figures A.3 and 


A.4@ illustrate the relationships. (Ref. Sj uses the 


parameters in table A.1l as the necessary elements of orbit. 


TABLE San! 
T - period (in minutes) w - argument of perogee 
i —.ngle of “inclination e - eccentricity 


Ae SSP longitude at perogee t - time at perogee 


The elements listed in table A.1 completely describe the 
theory of satellite motion. Perturbations are ignored. Im 
practice the motion of a satellite is very nearly described 


by these elements of orbit. 
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Figure A.1l Basic Ellipse 
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Figure A.2 Orbital Plane 
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Figure A.3 Elliptical Orbit Diagram 
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Figure A.4 Ground Tracks and Ascending Node 
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APPENDIX B 


RELATIVE MOTION 


The general theory of relative motion is developed using 
Simple frames of reference. An object in frame A moves with 
respect to that frame. Frame A moves with respect to the 
fixed frame B. These frames are considered to be coordinate 
systems or sets of reference axes. Figure B.1 depicts an 
object in the prime system. The X-Y axes represent the 
fixed system. In figure B.2 the prime system moves with 
respect to the fixed system a distance d. The object in the 
prime system moves a distance d’ within its system. 
Therefore the object moves a distance drersai. = d+d’ with a 
velocity Vreac = V+ V’. 

Further studies concerned with satellite tracking should 
involve relative motion. The fixed coordinate sytem is a 
three dimensional system with its origin at the center of 
the earth. The center of mass of the satellite is the 
Origin for the relative or prime coordinate system. The 
prime system orbits the fixed system. 

The target is detected by the satellite and its motion 
is therefore meassured with respect to the satellite’s prime 


coordinate system. These relative coordinates can be 
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related to the fixed coordinate system by a suitable matri 


transformation. 
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APPENDIX C 


ORIGINAL MODEL SIMULATION RESULTS 


This appendix contains the simulation for the original 
model. These graphs depict the states X and &. Figure 6.4 
presents the block diagram of the model simulated. The 
Simulations are run for 10 seconds and S00 seconds. This 
insures that a variety of suitable graphs are available for 
analysis. In order to accompany test observer behavior 
several different gains and time increments are used in the 
Simulations. Specific combinations of gain and time 
increment produce excellent results (as explained in 
chapter 6). 

The observer model’s ability to track a constant is 
recorded on eight graphs. The first two graphs (pages 112- 
113) depict state X for constant angular velocities and 
different gain values. The next two graphs (pages 114-115) 
depict state $ for the same gains mentioned above. A gain 
Ck) of ten is clearly beneficial for both X and &, since 
then Atk = O.1. 

The following two graphs (pages 116-117) show two 
different time increments for state X, constant gain and 
constant angular velocities. The graphs on pages 118 and 


119 display the same two time increments for state &. As 
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expected the smaller of the two time increments (0.01 
seconds) gives a slightly better performance. 

The last six graphs exhibit the model’s time response to 
Sinusoidal angular velocities. These are the results of the 
test used to determine the model’s ability to track time 
varying functions. The first three of these graphs (pages 
120-122) illustrate a target angular velocity of 
-OOl1lscos(Dt) and a satellite angular velocity of 
»YVOOS-cos(Dt). The variable D is varied in order to test 
observer response to different rates of change in angular 
VeloOecrtLes. ithe final graphs “pages 123-1235) depict X and & 
Vien Aueehaot angular velocity of -OOS*cos(Dt) and a 
satellite angular velocity of .00045-cos(Dt). The variabie 
Dis 0.01 in this simulation. This resultS in Sinusoids of 
smaller amplitudes and longer periods. The results of 
tracking the sine wave are very satisfactory. It 1S noted 


that the relative error for both states 1S very small. 
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APPENDIX D 


DECOUPLED M@DELS SUNULAT IGNBRESUETS 


These twenty graphs resulted from the decoupled model 
Simulations (Cas described in chapter 6). Figure 6.83 
presents a block diagram of the decoupled model. The first 
five graphs (pages 128-132) illustrate the time response tor 
state X. The time increment,.it, is held constant at 0.01 
seconds. Several different gains were used on each graph. 
The state X was followed closely by x. Different scales and 
Simulation times are plotted to produce & variety of 
Suitable graphs. A gain of ten consistently produces very 
good results in all five graphs. 

The next group of five computer graphs (pages 133-137) 
exhibit state ©. As above, the time increment, ‘t, is kept 
at 0.01 while the gain (k) is varied. The time response for 
the lower gains in particular are very satisfactory. The 
state $ is closely approximated by &. 

The third set of five graphs (pages 138-142) represent 
state X. This group of computer graphs holds the gain 
constant at one and allows the time increment to vary. As 
expected, smaller time increments (0.1 and 0.01) result in 
Superior tracking capabilities. 

The last set of five graphs (pages 143-147) show © and 


pee 


€. The gain is kept at a constant value equal to one and 
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the time increment is varied. As before, the smaller time 
increments produce the best results. 
In all of the above examples, the approximated state ON 


or $) closely tracks the actual state (X or 8). These 


Simulation results are quite satisfactory. 
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